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In Part I (in tliis journal) we argued that the structure function F2^{x,Q^) in deep inelastic 
ep scattering, regarded as a cross section for virtual j*p scattering, has a saturated Froissart- 
bounded form, evident already at HERA energies, behaving as ln^(l/a;) at small x, and that it 
can be extrapolated reliably to small x using the natural variable ln(l/a:). We showed further that 
the dominant structure functions F^^'"^ and FO'^^'^'' in charged- and neutral-current neutrino-nucleon 
scattering, respectively, can be related directly to Fi^^ up to the input of a small non-singlet Froissart- 
bounded function Tg. Here, we extend those results, and use them to evaluate ultra-high energy 
(UHE) cross sections for neutrino scattering on an isoscalar nucleon, A'^ = (n + p) /2, up to neutrino 
energies ~ 10^^-10^'' GeV where there are now limits on neutrino fluxes. The contribution 
of F2^'^'' is strongly dominant in the u-N cross sections, but we also include the contributions 
from the remaining structure functions xF!^^'^\x,Q'^) and Fj^^^\x,Q^), evaluated from F^^ with 
NLO corrections. We estimate that our calculated charged current and neutral current neutrino 
cross sections are accurate to ~2% at the highest energies considered (iJ^ = lO^'^ GeV), with the 
major uncertainty coming from the errors in the parameters that were needed to fit F2^ {x,Q'^) as 
measured at HERA. We compare our results to recently published neutrino cross sections derived 
from NLO parton distribution functions, which are very much larger than our cross sections at 
high energies because of the use of power-law extrapolations of quark distributions to small x. 
We conclude that our calculation of the UHE v-N cross sections is the best one can make, given 
the existing experimental deep inelastic scattering data. Further, the strong interaction Froissart 
bound of ln'^(l/a;) translates into a weak cross section bound of In^ E^, for v — N scattering, where 
Ev is the laboratory neutrino energy. Thus, the energy dependence oi v — N total cross section 
measurements have important implications for understanding hadronic interactions at enormous 
cms (center-of-mass) hadron- nucleon energies not otherwise accessible. 

PACS numbers: 13.15.-l-g, 13.60.Hb, 95.55.Vj, 96.50.S- 



I. INTRODUCTION 

Early in the development of perturbative QCD (pQCD), thepotential for dramatie growth of nueleon strueture 
functions as the Bjorken variable x became small was recognized Perturbative analyses showed that the number of 
low energy gluons rises rapidly as x, their fraction of nucleon energy, decreases When the collision energy is high 
enough, large numbers of small x quarks are generated in the parton sea, with interaction energies and momentum 
transfers large enough to be treated perturbatively. As a result, the collision cross sections of lepton, photon and 
hadron collisions on hadrons are predicted to show strong growth at ultra-high energies. This enhanced the prospects 
for detecting UHE neutrinos of cosmic origin, launched a number of related studies and has fundamental 
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implications for the design of exp erimental cosmic neutrino searches such as the past searches (AMANDA Q , ANITA 
[1,11, FORTE [H, GLUE RICE [H [H), those searches presently underway (ICECUBE [H, Baikal [H, 
ANTARES [Hi, HiRes [ll^UGER [H]), and those experiments under development (ARA [H], ARIANNA [H) or 
proposed (JEM-EUSO [2l|, [22[). All of these rely on theoretical models for total neutrino-nucleon cross sections cr'''^''^ 
at very high energies. Some searches among those listed already place limits on the neutrino flux; others are being 
designed for discovery of neutrinos above 10^^ GeV [2l|, [22|. Our work is designed in part with this energy target in 
mind. 

Much of the program at HERA, the electron-proton collider at DESY, was devoted to deep inelastic scattering 
(DIS) ep scattering experiments that probed small x in the range 10^^ to 10~^, with virtualities = — (e — e')^ 
that ranged from 0.1 GeV^ to 3000 GeV^. The exhaustive work by the HI and ZEUS detector groups confirms that 
the structure function F2^{x, Q^) rises rapidly as x decreases with fixed Q^. 

We argued in Part I of this paper [2^ that the structure function i^^^(x, Q^) is essentially the total cross section 
for the scattering of an off-shell gauge boson on the nucleon, a strong interaction process up to the initial and final 
gauge boson-quark couplings so we chose a fit function that saturated the Froissart bound |24| - [26| on total hadronic 
cross sections, (t(s) cx In^ s as the Mandelstam variable s for the gauge boson, nucleon system becomes large, s — > oo. 
This is suggested by the remarkably successful description of hadron-hadron and photon-hadron total cross sections 
over many orders of magnitude with the same saturated functional form in the relevant Mandelstam variables [27| . 
i.e., In^ s. Moreover, the very high energy proton-proton and proton-air cross sections predicted using this form for 
the LHC [28I - I30I I and the Pierre Auger Observatory [3l[, respectivel y, a re confirmed by the new data j32l - [34l |. Our 
high-quality saturated Froissart bound type fit to the small x data j35l - l37l |. extrapolated to ultra-small x using its 
known functional form, forms the basis of our present calculations. 

We also showed that the quark distributions at small x could be derived from our fit to -F^^ with only the additional 
input of a relatively small valence-quark contribution U — Uy ^ 2dy and the small non-singlet combination of quarks, 
Tg = u + u + d + d— 2s — 2s, which gives the difference between the light- and strange-quark distributions. We used 
the results to investigate the supposed "wee parton" limit, in which all quark distributions are assumed to converge 
to a common distribution at small x and large Q^. This limit had been used earlier [H, to calculate UHE neutrino 

cross sections in terms of -F;^^, neglecting contributions from the b quark and the neutrino structure functions F'^'^^^ 
and F^^'^^] the last vanishes in the wee limit for n/ = 4. We showed that the wee limit did not actually exist, but, 
surprisingly, that the deviation of the correct F!^^^^ from its wee parton limit was quite small, lending support to the 
approximate results in (sBI . l36j . 

Here we apply our results from Part I to a complete calculation of UHE v-N cross sections. We show first that the 
dominant structure functions F2^'^^ and FQ'^^'^^ in CC and NC neutrino scattering are equal to numerical multiples 
of up to small non-singlet corrections which we calculate to next-to-leading order (NLO) in the strong coupling 
as- We include the contributions of the b-quark and the subdominant structure functions F!^^'^^ and F'[''^\ F(f^'^\ 
calculated here including NLO corrections. We also extend the energy range of earlier calculations up to Ey = 10^^ 
GeV, the highest reach of the experimental search for UHE cosmic neutrinos [13 [m. Finally, we discuss how UHE 
neutrino measurements can give important information on hadronic interactions at energies far above those that are 
otherwise accessible. 

In the next Section we review the quark-parton construction of DIS e-p and v-N structure functions and the 
subsequent evaluation of the v-N charged and neutral current cross sections. We summarize the relations between 
ep and v-N structure functions from Part I, and apply the results to the representation of the small x charged and 
neutral current v-N structure functions in terms of F^^. We summarize the results of our very accurate fit to the 
HERA data on F^^(a;, Q^) in Sec. IIII Al and use the results as extrapolated to ultra-small x to make predictions for 
the UHE cross sections in Sec. IIII Bl 

In Section lill CI we show that our v — N cross sections are bounded by In^ F,^, where F^ is the laboratory neutrino 
(anti-neutrino) energy, a direct consequence of the hadronic scattering of the off-mass shell gauge boson on N being 
bounded by ln^(l/a;). In Sec. IIII Dl we discuss the implications of these neutrino results for UHE hadronic interactions, 
showing that an average hadronic cms energy W = a/I of ^ 70, 000 TeV is achieved in neutrino interactions with E^ 
of 10^^ GeV. We summarize and draw conclusions in Sec. IIVI 

II. ULTRA-HIGH ENERGY v-N CROSS SECTIONS 

A. Differential cross sections 

Expressions for the general quark parton charged current (CC) and neutral current (NC) v-N cross sections are 
given in many references; see for example [gI. |38| - |40| . We display them here to keep our presentation self-contained, 
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using the notation of [39[ . The double differential inclusive charged current cross section for neutrino or antineutrino 
scattering on an isoscalar nucleon target N = {n + p)/2, i^i + N—^i + X or Dg + N^i + X,£ = e, /i, r, is 



47r 
1 



X ■ 



X 



±xF: 



\ 2nixEu 



(1) 



The upper signs are for v and the lower, for v. In this expression, x = Q^/2p • g is the Bjorken scaling variable where p 
is the nucleon 4 momentum, and = —Q^ where q = v — t\s the momentum of the virtual W boson which interacts 
with the nucleon, that is, the momentum transferred from the leptons in the scattering. The second independent 
scaling variable is the fractional energy transfer to the nucleon, y = {E^ — Ei)/E^\ clearly, < y < 1. In covariant 
form, y — p ■ q/p ■ v = /2mxE^, as it appears in Eq. ([T]). 

The direct channel v~N Mandelstam variable s is s = (i^ = 2'mE^, m the nucleon mass, where we neglect 

m? relative to 2mE^, so y = Q^/xs. Also, the direct channel Mandelstam variable s for the W*N scattering is 



s = (q+pr ^2p-q-Q^ 



Q /x. The structure functions Fi, F2, F^, and F^ ^ F2 — 2xFi are functions of the 



Bjorken variable x and the virtuality of the gauge boson. 

The functions Fi = Fi, x^^F2, F3 are given in terms of quark distributions by expressions Fi^ of the LO form 
convoluted with QCD correction terms |4l| - |43| . schematically 



F, 



1 + 



ZTT 



(2) 



Here / is the unit operator, Ciq and Cig are coefficient functions from the operator product expansion, known in 
low order. The symbol indicates convolution of Ciq, Cig with Fi^o and the giuon distribution g = g(x,Q^) = 
x~^G{x, Q^), where the convolution of operators A and B on x space is defined as 



f dz f dz 

A®B= —A{x/z)B{z)= —A{z)B{x/z). 

The uncorrected structure functions Fi q are sums of quark distributions 



(3) 



i j qj (a;, Q ) , 



(4) 



with the constants Cij being the weights with which the different quarks appear. The parton distribution functions 
and coefficient functions are all to be calculated at the relevant order in the strong coupling as- More detailed 
expressions are given below. 

To express the charged current cross section in terms of PDFs, we write the structure hmctions 
Fi q{x,Q'^), F2^q{x,Q'^) and xF^q{x,Q^) for neutrino scattering on an isoscalar nucleon N = {p + n)/2 in terms 
of the quark distributions u = u{x, Q^), d = d{x, Q^), s = s{x, Q^), . . ., as follows: 
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In these expressions, we have referred all distributions to the proton so that c?„ — ^ Up 



u, Un 



d, etc. We 



will take s = s, c = c, b = b since these quarks are produced only in pairs through gluon splitting. It is then clear 
that F2Q = F2Q and that F^q and F^q would have opposite signs except for the presence of the valence distributions 
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: d ~ d. Valence-quark effects are unimportant at small x, the region of primary interest here, so 

The t quark, taken as active as a parton only for > m^, does not contribute significantly to the high-energy 
neutrino cross sections because of the effective cutoff in the integration over for ^ -^w^ ^'^"^ contribution 
will be taken to be zero in the present calculations. The b quark, active for > ml, contributes to the CC cross 
section only when the threshold condition s k, Q"^ /x > for the process b ^ t is satisfied. 

The uncorrected longitudinal structure function Fj^ q — F2Q — 2xF^q is equal to zero in leading order in the strong 
coupling, but the corrected function F^j^ gains terms of orders as /(27r) and higher through the convolutions of the 
different coefficient functions for F^ and a;~^F^ with the quark and gluon PDFs [isj . 

The corresponding expressions for the neutrino neutral current structure functions and cross sections are 
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= {d + u + d + u) {Ll + Rl + Ll + Rf) /4 

+ {s + b+~s + b){Ll + Rj) /2 + {c+ -c) {Ll + Rl) /2, 
= ±{d + u-d~u){Ll-Rl + Ll-Rl) /A, 



(11) 

(12) 
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where FO 



denotes the neutral current structure function. The NC chiral coefficients are defined as L,, 
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2 oiii lyj/v , -ufl — 1 -|- ^ sin^ ) Ru — — ^ sin'^ and i?^ = ^ sni c^y, 
calculations. The double differential neutral current cross section is 



I sin^ ' 



0w, and the value sin 6'v^=0.231 l40| is used in all 
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B. The dominant structure functions ^2"^''' at NLO 

For reference, we recall that the uncorrected structure function x^^F^q electron-proton DIS is defined in terms 
of quark distributions qi — qi{x, Q^) of a given order as 



(15) 



where the e,; are the quark electromagnetic charges in units of e, and the sum is over the active quarks and antiquarks. 
This form gives the complete Fj^^ at LO. 

Beyond leading order, F^^ is given by (4l| - |43| 



x-'f: 



1 p7P 



(16) 



where 1 is the unit operator and (8) indicates convolution of the coefficient functions C2q and C2g with x ^F^q and 
the gluon distribution g, as defined in Eq. 

The coefficient functions C2q and C'2g depend on the renomalization scheme used in perturbative calculations; we 
assume the use of the standard MS scheme in which, at NLO [ij . 
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The expression in Eq. p6|) is usually used to determine F^'' from the individual quark and gluon distributions 
found in fits to the DIS data, as, for example, in [13, However, the relation can also be inverted to determine F^q 
directly at a given order in as in terms of F^'' and a given gluon PDF. This is discussed in detail in the Appendix to 
[2^. We will assume that F^q is known. 

As shown in Part I [1^, the neutrino structure functions F^q'^'^ can be expressed in terms of F^q, the valence 
distributions ~ u — u and = d — d, and a set of non-singlet combinations of quark distributions T^2_i. Since 

2 



above the 6-quark excitation threshold 



valence effects are unimportant at small x, we will use the reasonable approximation d„ (l/2)C/i, = U . Then for Q 

F20 = x{u + d+2s + 2b + u + d+2c) (19) 
= YI^2lf-^(5r8-5ri5 + 3T24)-iJ 



-F20 - (57^8 - 5ri5 + 3T24) - —U, (20) 



where we have ignored the very small differences between the u and d distributions at small x. The Ts are given in 
terms of quark distributions by [43| 

Ts = x{u + u + d + d-2s~2s), (21) 
Ti5 = x{u + u + d + d + s + s-3c~3c), (22) 
T24 = x{u + u + d + d + s + s + c + c-4b~4:b). (23) 

Wc showed in Part I how T15 and T24 can be constructed given F2Q and the input only of U, taken from CTEQ5, and 
Tgfcc, m^), which we took as the extrapolation to small x of the result given in NNLO by the CTIO quark distributions 



47[, using an extension quadratic in lii^^l/x) as required by the form of our Froissart bounded fit to F^^. Since 



F2Q can be determined, and the evolution of the Ts is minimal and known in NLO |23j . we will take the Ts as known 



and will use them instead of the individual quark distribution in the following discussion. 

We obtain the observable neutrino structure function Tj"''^'' by applying the QCD corrections from the operator 
product expansion [Hl-lillii] to F2^^\ 

x-'f!:^"^ 
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where the coefficient functions C2q and C2g are given to NLO in Eq. ([T7)) and Eq. ([T5| . and the C2,i are the coefficients 
of the quark distributions in Eq. ([T^ . This transformation converts the T^q^ term in Eq. ([^0)1 to Tjj^ as in Eq. (|16p . 
The gluon term in Eq. p4)) is absorbed in the process, and only C2q acts on the T,- in Eq. ((20| with T,- Tj = 

x [1 + (as/2TT)C2q] €5 (x~^Tj). This transformation of the Ts is discussed in Part I. The result for F2^^'' for different 
numbers of active quarks is 

^U', nf = 3, Ql<Q^<ml, (25) 

^It[,-^U', nf = A, ml<Q^<ml (26) 

+ ^Ti,-^TU~^U', nj = 5, ml<Q\ (27) 

In the wee parton picture, the quark distributions all converge toward a common distribution at small x where 
the valence quark contributions are negligible and sea quark distributions are all equal (equipartition of flavors), 
xqi{x, Q^) — >■ xq{x, Q^) for all i. The Ts then vanish individually, and 

F^^'^ « If^, nj = 3, Ql<Q'< ml) (28) 

Fl^^'-'^^Fr, nf=A, ml<Q'<ml (29) 

F2^'^ « ^F^, n/ = 5, ml<Q\ (30) 

The coefficients in these expressions are just the ratios {J2i C2,i) / {J2i ) of the sums of the squares of the weak and 
electromagnetic charges, summed over the active quarks. 

The wee parton picture was used in Refs. [ssl. fSal with n/ = 4 to predict UHE neutrino cross sections in terms of 
Tj^, neglecting the contribution of the b quark and the small terms that involve the structure function F^'''^\ The wee 
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parton limit actually does not exist: as wc showed in Part I, the quark PDFs for different flavors diverge from each 
other rather than converging at small x, with the ratios determined by the x dependence of f^^^. However, we showed 
also that the combination of T's in Eq. (|20|) is quite small relative to F^J' because of cancellations, and decreases the 
final result by ^5% to 3.6% as x decreases from 10~^ to 10~^^ at = 100 GeV^, and by 1.5% to 1% over the same 
X range for = 10, 000 GeV^. The smallness of the combinations of T' and U' terms leads to effectively the same 
results as obtained in the wee parton limit. 

The manipulations which connect F^^'^'' to F^'^, culminating in Eqs. (pSjl - ipT]) . are discussed in detail in [2^. The 
final results for the leading terms are correct formally to all orders in as , and F2^ is known to experimental accuracy 
from our fit to the HERA data. The theoretical uncertainties in the final result arise through the non-singlet functions 
U', Tg, T^'g and T24, which we have only treated to NLO in their evolution and under the transformation in Eq. (|24p . 

In our calculations, the valence distribution U was taken from CTEQ5 (49| and the (small) function Tg was taken 
from the result from the CTIO analysis of the HERA and other data [4y|, extrapolated to small x using the Froissart 
bound form of the fit function required by consistency with the form of F2^ . They were used as an input to get the 
singlet distribution Fg — + 9i) from F2Q, known to NLO, at = m^. This determines 115(2;, m^) = Fs{x, ni^). 

As shown in Part I, the changes in Tg and T15 induced by QCD evolution to = to^ are minimal and can be 
calculated analytically to sufficient accuracy for our purposes, and T2a{x, ml) then determined from F;7cf , Tg, and T15; 
the evolution of T24 is again known analytically. There is a further hidden uncertainty in that the transformation 
from F2^ to the uncorrected distribution F2Q needed in the determination of Fs depends on the gluon distribution t/, 
which we again took from the CTIO analysis as extrapolated to small x. We estimate the overall uncertainties in the 
small non-singlet T terms in Eqs. (|25p - ([27)) to be significantly less than the < 30% difference between F2^ and our 

derived F2Q , and to lead to at most a 1-2% uncertainty in the final results for F2'''^\ 
The analog of Eq. ^ for i^Oa^^'^for > ml and n/ = 5 is 

FO;^'^ = Y2 {^"^' 1^^') [2 + + 2 [Li + Rl)] 

[4 {Ll + Rl) - {Ll + Rl)] (5T^ - 5T[, + ST^) ■ (31) 

The contributions associated with the T's are negative and decrease the final results for FOj^"^^ by ^5.9% to 4.3% 
as X decreases from 10^^ to 10^^^ at = 100 GcV^. and by 1.8% to 1.1% over the same x range for Q'^ = 10,000 
GeV2. 

Neglect of the T' terms, which would vanish for all quark PDFs equal, again gives an effective wee parton limit, 
valid to a few percent even though the limit does not exist at the level of the individual quark distributions. For 
completeness we write out these approximate relations for all nf. 

FO^^'^ « ^Fr [Ll + RI + 2 {Ll + Rl)] , Ql < < ml (32) 

FO;^'^ - ^^^2^^ [Ll + RI + LI + Rl] , ml<Q'< ml (33) 

FO^^'^ - ^F;'^ [3 {Ll + Rl)+2 {Ll + Rl)] , ml < Q\ (34) 

In our previous work [111 and [111, the electron-proton structure function F2^{x,Q'^) was provided by fitting the 
early ZEUS data [111 and later, the combined ZEUS and HI analysis of the full HERA data set [sij. We will use the 
latter fit in the present work, which goes beyond the earlier papers not only by establishing the link between .F^^ and 

F^'"^^ in general, but also by calculating the complete v~N cross section including the non-singlet corrections and the 
structure functions F'^^'^^ and F^'"^-* in NLO, including the b-quark contribution which appears through T24. 

C. The sub-dominant structure functions F^^"^ and 7^^^"^ at NLO 

Although the F!^'^^^ structure function is dominant in v-N and D~N UHE scattering as will be discussed later, the 
contribution of the structure function F!^^"^ through the b quark is significant; we wish to estimate it at NLO, and 
further, to include the NLO contribution of F^^''\ which is zero at LO in a^. We will concentrate on the region 
> ml which contributes all but a small fraction of the cross sections, as we will see in Sec. IIIII 
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We start by re-expressing xF^^"\ Eq. ([7]), in terms of -F^o''^ ^^"^ ^^"^ 
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11, 



n/ = 3, 



9 pp(p) 



20 
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STis + 3r24) + — c/. 



(35) 
(36) 
(37) 



The expression for Uf 



5 has the same structure as that for F^^^^ in Eq. (|20p , but with a coefficient for the T term 
which is larger relative to that for the i^2cf term by a factor 25/3, and a valence term of the opposite sign. As a 
result, the combination of Ts and J7, previously quite small relative to i^^J', is now significant and should be taken 
into account. If it is not, we get as the effective wee parton limit 



30 



(9/11)^2^^ 



u{y) 



(38) 



a relation useful for making rough estimates of K 



To obtain the physical structure function K 



30 



we must convolute F^'^ with the QCD coefficient function C^q, 



^p(p) _ _pp(p) 



3,0 



277 • 
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pP(p) 
3,0 



(39) 



This is given in NLO by 



C3,(z) = C2,(z)--(l + z), 



(40) 



with C2q{z) the function in Eg. p6|) : see Refs. [4ll443| . The final structure function xF^^'^^ depends only on a 
convolution involving the quark fields; the gluon does not enter and no assumption about g is necessary even though 

o"^ has a mixed singlet, non-singlet structure for = 6 nonzero Since the gluon does not appear and C^q ^ C2q, 
the transformation does not convert F^n to F2^ as in the calculations above. However, the convolution integrals are 



20 



readily evaluated at small x using the methods of Part I [23[ once xFj^q'^'' is known; see the Appendix of Part I. 

The contribution of the resulting xF^''"'' to the cross sections is described in Sec. IIIIl where the cross sections are 
evaluated and the various parts of the final results compared. We note here that at small x, where the valence term 
U is negligible, xF^{x,Q^) = —xF^{x,Q^). As a result, a^^ « a^i^ for large in both LO and NLO despite the 
presence of the ± signs in Eq. ([T]). 

U and 



The NC structure function FO'^q''\ Eq. depends only on the valence quark distributions 



dv = d 
there. 



d ~ U/2, so vanishes at small x where u ^ u and d d, and docs not contribute to the NC cross section 

pp(p) 



The structure function Fj^ '{x, Q ), which is zero in LO, is given in NLO [43| by 



x-'F^^^'\x,> 



) = ^CLq®ix- 



20 



7^2nf CLg ' 



CLgiz)^2zil~z). 



(41) 
(42) 



The integrals can be calculated analytically for small x following the methods outlined in the Appendix to [23[. It 
is sufficient for our purposes to approximate -F^j'o'^'' '^^ this calculation by {45/11)F2q for nf = 5, as in Eq. We 
include F^^'^-' in the cross section calculations in Sec. IIIIl 

A similar result holds for the NC structure function F(f^^\x,Q^), with F!^^^^ replaced in Eq. (PT|) by FQ 

{%/ll)F^'' - 



[Ll + Rl + {i/2){Ll + Rl)l¥.^.m- 



-|p(p) 

-"20 



D. Total cross sections at UHE 



Because the UHE total cross sections for neutrino-isonucleon scattering are so important for neutrino observatories 
and for testing for new physics beyond the realm of laboratory experiments, we lay out the details of our calculations 
and comment on several of their features. 



8 



The inclusive double diiTerential cross sections in Eqs.(lT|) and ([T4|) for the inclusive CC and NC cross sections for 
the processes I'l + N ^ £ + X and + N ^ i^f + X , £ = e, ^, t, depend explicitly on the invariants s = {v + p)^, 
= —q^, q ~ V — I, V — v' , and p ■ q through the scaling variables x — Q^/2p ■ q and y = p ■ q/p ■ v = /xs. Since 
p\ = {p — q)^ > m^, m the nucleon mass, — + 2p ■ q > and a; < 1. Also, if E'l, ^ to in the nucleon rest frame, 
then s = 2mE^, and y = {E^ — E'^,)/ in that frame. Thus y is the fraction of the neutrino energy transferred from 
the incoming neutrino to the outgoing lepton, with < y < 1. As a result, the relation = 2mE^xy requires that 
< 2raEy — s. The structure functions p'^''^^ are functions of x and . 

We incorporate these conditions in the integrations over and x. We limit ourselves to a minimum value of that 
is consistent with the application of perturbative QCD and at the same time avoids possible problems with numerical 
instabilities as x or goes to zero. We therefore take and x in the ranges s = 2mE^ > > Q% i„ — 1 GeV^ 
and / {2mE^) < x < 1, integrating first over x. As emphasized in early calculations (see Ref. jS^I), the vector 
boson (V) propagator factor {My /{My + Q^))^ cuts off the integrand for >> My-, effectively selecting a range of 
small X reaching somewhat below x ~ My / (2mEi,) which makes the only substantial contributions to the total cross 
section. For the range of neutrino energies we consider in this work. 10^ GeV < E^ < 10^^ GeV, this means that we 
must accurately probe x values in the range 0.001 < x < 10"^'^. 

With the above discussion in mind, we express the total cross sections as 



'CC 



dx 



47r . n2 -1 



dxdQ^ 



Mir 



dx 



f:^ + xF^ + {f!^^xF!^){i-:^\ ~[^\fi 



Q2/(2ni-B„) ^ 

2X2 /g2x2 



XS 



xs 



(43) 



The corresponding neutral current total cross section, a'^Q{E^) is obtained by the replacement of Mw by Mz and 
the replacement of F^ , xF^ and F^ by FO^ 



xFO^ and FOl. 



Written in this form, the integration over x shows clearly that a simple power law behavior of F2^'^'' as a function 
of X (or of ln(l/a;)) will produce the same power of x (an added power of ln(l/a;)) in the result. This is helpful for 
assessing the UHE behavior of the total cross section that follows from a given model of the structure functions, as 
pointed out in [s^]. In particular, the ly-N cross section calculated to lowest order in Gp, as done here, will rise 
asymptotically as ln^(l/a;) cx Ib!^ E^, for E^/Q'^ ^ 1 for our Froissart-bounded extrapolations of F^*-"-*. However, the 
authors of [l^ incorrectly state that the work in [s^ and [s^ claims that the v-N cross section rises as In^ E^, in the 
large E^, limit. This confusion between the ly-N and gauge boson- iV cross sections is straightened out in (53 |. 

III. PREDICTIONS FOR ULTRA-HIGH ENERGY NEUTRINO CROSS SECTIONS 
A. Extrapolation of F^'' to ultra-small x 



The expressions in the preceding sections together with the methods developed in Part I allow us to calculate the 
main contributions to the v-N cross sections directly in terms of F^''. However, as we emphasized in the Introduction, 
for the energies Ei, ~ 10^^ GeV of interest for UHE neutrino cross sections, we must know F^^(x, Q^) at values of x 
6-7 orders of magnitude below the lowest values ~ 10^^ encountered at HERA, or at values of the natural variable 
V = ln(l/a;) a factor of ^ 2-2.5 larger than at HERA. 

Given these conditions, it is essential that the form used to extrapolate F^^(a;,(5^) be consistent both with the 
asymptotic limiting behavior expected theoretically, and with the present HERA data. We can, in fact, only expect 
the extrapolation to be unambiguous and robust if an accurate fit to the data indicates that the measured structure 
function is already consistent with the expected asymptotic form. This is the situation we encounter here; the resulting 
extrapolation is the best one can do. 

In Part I of this work [H], we reviewed and strengthened the case that F2^{x, Q"^) saturates the Froissart bound. We 
adopt the global fit function used in (5^, [s^ and 37|, which insures that the Froissart-bounded behavior, ln^(l/x). 
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dominates at small x. It takes the form 
where 



Fp 



Xp 
Xp 1 



^(Q')ln 



Xp 1 



X 1 — Xp 



X 1 — Xp 



(44) 



hi In 



(45) 



The expression in Eq. (|44p is equivalent to a quadratic in the logarithm of the square of the 7*p center of energy 
s — ~ Q'^{1~ x)/x with the dependence rearranged and made more flexible. A quadratic in Ins, s the square 
of the hadronic energy in the center of mass system, is known to describe hadronic and real 7-nucleon scattering to 
the highest energies studied [l^, l32| - [33 | . It is therefore natural to use it to describe what we regard here as a hadronic 
process, 7*73 scattering, while allowing through the coefficients for the extra dependence introduced by the mass 
of the virtual photon. 

The procedure used in fitting the combined HERA data is described in [3^[33|. We fixed the parameter xp at the 
value 0.11, and and used the functional form in Eq. (|44)) to fit the combined HERA the data for x < 0.1; there are 
very few HERA data for larger x. Fp, the common value of i^^^ at xp, and the other 6 parameters are shown with 
their errors in Table HI To extend the results for F^^ to x > 0.11, wc used the extrapolation method discussed in [37j . 

The renormalized minimized value [56j , the number of degrees of freedom and the rcnormalizcd Pcr degree 
of freedom are 391.4, 335 and 1.17, respectivel y, f or our analytic global fit [s^ to the combined ZEUS and HI results 
[Slj . We provide a fuller discussion in Part I jl^. Since the fit depends only linearly on the parameters, the errors 



Parameter 


ao 




a2 


bo 


61 


&2 


Fp 


Value 


-8.471 X 10"^ 


4.190 X 10-2 


-3.976 X 10-3 


1.292 X 10-2 


2.473 X 10-* 


1.642 X 10-3 


0.413 


Error 


±2.62 X 10"^ 


±1.56 X IQ-^ 


±2.13 X 10-" 


±3.62 X IQ-* 


±2.46 X IQ-* 


±5.52 X IQ-^ 


±0.003 



TABLE I: Parameter values and errors for the global Froissart bound model fit function. 



propagate linearly in our Froissart bound type model, which means that the percentage errors in the results of our 
ultra-small x extrapolations should be comparable to those of the parameters themselves. 



B. Results for ultra-high energy neutrino cross sections 

Given our fit for i^^^(x, Q^), it is straightforward to calculate the expected neutrino cross sections by integrating 
Eq. ([1]) or Ec^. over x and Q^, as in Eq. P5|) . We consider the case of five active quark fiavors. Because of the 
effective cutoff imposed by the gauge-boson propagator factor l/((5^ ± My)'^ in the integrand, the top quark, 
which becomes active as a parton only for >> w (173)^ GeV^, docs not contribute significantly. However, as 
we note later, the real CC process + & — s- 1 is allowed over essentially all of the important region in x, space. 

The simplest calculation, accurate to a few percent at ultra- high neutrino energies E^, , uses the effective wee parton 
approximation to express the structure functions in Eqs. ([1]) and (jl4p directly in terms of F2^. However, it is not 
much more difficult using the methods of Part I to use the full expressions for those quantities. The additional terms 
are generally small, and it is presumably sufficient to treat them to NLO, as we have done. The general accuracy of 
the wee approximation is shown in Table HIl where we compare the two calculations. 

In Fig. [T] we show the charged current and neutral current neutrino cross sections for large E^, which follow from 
our global Froissart-bounded fit to F2^ using our complete results. We believe these give the best predictions that 
can be made using current experimental information. They do not depend on the large extrapolations of the starting 
PDFs used in standard analyses — typically with power dependence in x — over many orders of magnitude in 1/x; they 
depend only on the extension of our physically motivated Froissart-bounded fit to the HERA data on F2^ over a 
factor of « 2.5 in the natural variable v = ln {1/x). 

As the neutrino energy increases, the interaction probes more deeply into the small x region of the nucleon. We 
illustrate this effect in Fig. [51 where we show the differential distribution of contributions to the CC cross section with 
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Energy (GeV) 


10'' 


10* 


10^" 


10^2 


10'^ 


10^^ 


CrCC, exact, Cm^ 


7.40 X 10"=^^ 


4.89 X 10"^^ 


1.82 X 10"^^ 


4.69 X 10"^^ 


9.75 X 10-32 


1.76 X 10-3^ 


crcc, woo, cm^ 


7.54 X lO-^*" 


5.00 X 10"^^ 


1.84 X 10"^^ 


4.74 X 10"^^ 


9.84 X 10^32 


1.78 X 10-31 


CTNC, oxact, Cm^ 


3.04 X 10"^'' 


1.97 X 10"^^ 


7.58 X 10"^^ 


2.00 X 10-32 


4.20 X 10-32 


7.67 X 10-32 


2 

o"NC, woo, cm 


3.12 X 10-=*" 


2.02 X 10"^^ 


7.70 X 10"^^ 


2.02 X 10^32 


4.24 X 10^32 


7.73 X 10^32 



TABLE II: Comparison of the exact CC and NC cross sections with their wee parton approximations. 



1x10"^' 
5x10-^2 



B ix 10-32 



1x10-3 

5x10-3 

1q6 jq8 iQlO io'2 IQW lQl6 

£v(GeV) 

FIG. 1: Plots of the u-N cross sections, in cm2, vs. the laboratory neutrino energy, in GeV, calculated using the small x 
extrapolation of the global fit to the HERA data on F^^ix, Q^) given in Eq. (|44|l utilizing the relations between F!^'^"^ i^Oj'''^ 
and F2^ in Eqs. H27|) and (|31|l . with a NLO treatment of the small functions T/ and of the subdominant structure functions 
F^^''^ and F^'"'. The upper curve (red) is the CC cross section and lower curve (black) is the NC cross section. 




respect to a;, normalized to the total cross section. The dominant contributions to a'^^ {E^) march steadily downward 
in X as Ey increases. For example, at Ei, = 10^^ GeV, the region in x from 10^® to 10^* determines the value of 
^"""{E,)- ... 

Given this information, we see that our extrapolation of F^^ into the ultra-small x region is much less extreme than 
it might be, since the Froissart-likc form for i^^^(x, Q^) reduces to a quadratic form in w = ln(l/a;) for x small. The 
bulk of HERA data arc in the range from v ~ 2.3 to w w 10. It is evident from Fig. [2] that the most important region 
for the prediction of the neutrino cross section at Ey = 10^^ GeV is around x = 10^^^ or v = 25, an extrapolation 
of roughly a factor of 2.5 in this variable. The fact that the small x or large v form of -F^^'j ^ simple quadratic in 
V, is tightly constrained theoretically by the Froissart limit, and experimentally by the quality of the fit, limits the 
uncertainties introduced by the extrapolation. 

The values of the CC and NC cross sections and their fractional uncertainties calculated from the squared error 
matrix of the fit (including correlation errors) over the energy range from 10^ GeV up to 10^^ GeV are shown in 
Table IIIIl The integration errors are negligible, a part in 10^ or better. As is expected for a linear fit, calculated 
uncertainties of 1% to 2% are similar is size to the errors in the fit parameters themselves. As far as purely numerical 
uncertainties are concerned, these are the smallest one can obtain within the range of current efforts to estimate the 
UHE CC and NC neutrino-nucleon cross sections. 

The dependence of the I'-N CC cross section is shown in Fig. [31 where we plot the partial cross sections obtained 
by integrating Eq. ([T]) first over x, and then over with > Qq, for a selection of values of Qq from 1 GeV^ to 
s = 2mEy. The partial cross sections are normalized by the total CC cross section and plotted against Qq, the 
minimum value of included. The curves show the results for E^, = 10^, 10^, 10^", 10^^, 10^** and 10"^^ GeV. 
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10-11 10-' 10-^ 10-5 0.001 



X 

FIG. 2: The normalized distributions of contributions to the u-N total cross CC section of Fig. [T] with respect to x for cross 
sections at energies 10*, 10^", 10^^, 10^'* and 10^'' GeV, right to left. The integration over the variable x involves dx/x, so the 
In a; distribution is the natural measure. 



Energy (GeV) 


1Q« 


10^ 


10« 


10^ 


10i« 


10" 


10^2 


10" 


10" 


10" 


10" 


10^^ 


(Tcc (nb) 


0.740 


2.06 


4.89 


10.0 


18.2 


30.2 


46.9 


69.1 


97.5 


133 


176 


228 


Sacc/o'cc 


0.009 


0.010 


0.012 


0.014 


0.016 


0.017 


0.017 


0.014 


0.015 


0.015 


0.016 


0.017 


CTNC (nb) 


0.304 


0.817 


1.97 


4.12 


7.58 


12.7 


20.0 


29.6 


42.0 


57.6 


76.7 


99.6 


^CTNc/fNC 


0.010 


0.010 


0.013 


0.015 


0.016 


0.017 


0.017 


0.014 


0.015 


0.016 


0.016 


0.017 



TABLE III: Total cross sections in nb and the correlated fractional errors in the v-N CC and NC cross sections, as a function 
of neutrino energy in GeV, computed from the global saturated Froissart-bounded fit to the F^^{x, Q^) HERA data. 

Combining these curves with those in Fig. [5J we can now identify the major contributors to the total cross section 
in Eq. (|43|) . We note that ii y = / {2xmE,j) « 0, only the structure function contributes to the cross section, 
with a coefficient of 2. As an example, at E^, = 10^^ GeV, we see from Fig.[5]that the most likely value of a; is « 10~^, 
and from Fig.[3l that sa 90% of the cross section comes from smaller than about 4 x 10^ GeV^. Thus, about 90% 
of the time, y < 0.2; so the coefficient of Fj" > 1.6. A complete calculation shows that the average value of y over the 
distribution ranges from ^ 0.2 at E^, = 10^ GeV to ^ 0.08 at 10^^ GeV. Similarly, we estimate the typical coefficient 
of F^^''\ itself < {1/5)F^^''\ as < 0.4, for a contribution relative to that of F^"" of - 0.05. The coefficient of the NLO 
function whose evaluation requires a knowledge of the gluon distribution in addition to F^q, is ^ 0.04. Thus, 
we see that the overwhelmingly dominant contribution to the neutrino cross section is from F2 ■ 

As shown in Table HIl F^'-'^-' is given to good accuracy as a simple multiple of F^*'. Further, as shown in Part I and 
discussed in Sec. Ill Bl this relation holds to all orders in ag. The relations of the non-singlet corrections T' and the 
sub-dominant functions F3 and F^*-"-* to are known to NLO, sufficient accuracy for our purposes since their 
contributions to the cross section are either intrinsically small, or are suppressed by the factors estimated above; all 
are included in our final cross section calculations. The CC v-N cross section is therefore determined to high accuracy 
by , a quantity obtained from the fit to experimental data. Simlar arguments hold for the NC cross sections. 

We also see from Fig. [3] that at least 90% of the cross section is captured in all cases for a lower cutoff in the 
Q2 integration of Ql = 500 GeV^, and, similarly, for Ql = 1000 GeV^ for E^ > 10^ GeV. The final state hadronic 
invariant mass ^ /x ^ is therefore greater than rnj over all of the important region oi x, space in 
the ultra-high energy CC and NC cross section integrands, e.g., for x < 0.017 for = 500 GeV^. The process 
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+ 6 — >■ t is therefore allowed and the &-quark contribution is essential for the calculation of the ultra-high energy 
CC and NC cross sections. This contribution was omitted in early work. 




1 10 100 1000 10* 10' 

Ql 

FIG. 3: The normalized distributions of contributions to the u-N total cross CC section of Fig.[3]with respect to Qq, in GeV^, 
for cross sections at energies 10^, 10*, 10^", 10^^ 10^* and 10^*^ GeV, bottom to top. After integration over x, which starts 
at /s, the integration over the variable starts at the minimum value Qq consistent with the range of validity of the 
determination of F2 {x,Q^). The plot shows that the calculated value of the UHE cross section is very insensitive to Qq, the 
lowest value of the integration range over Q^ . 

In Fig. |4] we compare our UHE cross sections from Fig. [l]with those of Cooper-Sarkar, Mertsch, and Sarkar [3^, 
who used the HERA-based PDF set HERAPDF1.5, and included the b-quark contribution in their computations. 
Their quoted error estimates are in the 2%-4% range, comparable to ours, when they exclude those PDF sets which 
lead to unacccptably steep rise in the cross section or allow negative values of the gluon PDF at small x and small 
Q\ 

The Froissart bound based and PDF based calculations agree very well for E,^ k, 10^-10^ GeV, where essentially 
the entire neutrino cross section arises from regions in x (Fig. [2|) and (Fig- El) corresponding to the x, region 
of the HERA data where the 5-quark is reasonably above its excitation threshold and valence-quark contributions to 
are small. 

At Eu = 10^^ GeV and 5 x 10^^ GeV, the highest energies reported in Ref. [s^, their cross sections are a factor two 
larger than ours and extrapolate to more than an order of magnitude larger than ours at E^ = 10^^ GeV. This large 
difference results from their use of a power-law extrapolation in x of the HERAPDF1.5 parton distributions they use, 
whereas our partons are bounded by the saturated Froissart-bound fit to F^^(a;,(5^), and grow only as ln^(l/a;). A 
measurement of the CC cross section for E^, > few x 10^^ GeV could provide a crucial test of our results. 

The application of solutions to the DGLAP equations for quark PDFs to evaluation of UHE ly-N CC and NC cross 
sections has a long history. The application reported by Ghandi, Quigg, Reno and Sarcovic (GQRS) @, based on 
d, u, s and c quark PDFs from 1998 CTEQ4 analysis of the early ZEUS small x data, was a standard for many years, 
and still remains a point of comparison. The recent results of Connolly. Thorne and Waters (CTW) (sj, and those 
of Cooper-Sarkar, Mertsch and Sarkar (CSMS) [1^ shown graphically in Fig. 01 include the 6-quark contribution to 
both CC and NC scattering, and are based on updated PDFs derived from newer and larger data sets. To provide 
points of comparison, we re-tabulate our results together with those just mentioned in Tables IIVI and IVl The tables 
cover the cross section up to the E^, = 10^^ GeV values published in GQRS and CTW. Since CSMS quotes cross 
section values only up to = 5 x 10^^ GeV, we enter them at E^ = 10^^ and note them with an asterisk. As already 
noted, CSMS use PDFs from a fit to the combined HERA results, so the data used and the inclusion of the 6-quark 
make their work the most natural to compare to ours. CTW use the MSTW2008 NLO PDFs, which bases its small 
X information on ZEUS data and provides PDF grids down to x = 10~^; for smaller x, CTW extrapolate the quark 
PDFs with a form a + blii{l/x), a and b constants, rather than the power-law form that follows from the extension of 
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5x10-3- 

^ 1x10-3- 
I 5x10-33 

b 

1x10-33 

5x10-34 



lO*" 10' 10** 10' 10'" 10" 10'2 

£v(GeV) 

FIG. 4: Plots of u-N cross sections, in cm^, vs. tiie laboratory neutrino energy, in GeV. Our CC cross section is the upper 
solid (red) curve and our NC cross section is the lower solid (black) curve; the CC cross section of Ref. [s^ is the upper dashed 
(red) curve and their NC cross section is the lower dashed (black) curve. The energy range is that reported in Ref. [s^. All 
cross section calculations include the b-quark. 

the initial MSTW quark parametrizations to small x. Use of the latter would lead to much larger cross sections than 
CTW obtained. 

All these PDF calculations use extrapolations of individual quark distributions to values of x well below the region 
of the HERA data using somewhat arbitrary assumptions about their x dependence. Our extrapolation, in contrast, 
is of the Froissart-bounded form of the full structure function - that is, the virtual-boson, hadron scattering 
cross section - which provides excellent fits to the j*p HERA data, and the jp and hadronic data to the highest 
energies studied. We regard this extrapolation, by a factor of ~ 2.5 in the natural variable v = ln(l/.T), as far more 
reliable theoretically and numerically. We conclude that the cross sections predicted by the PDF-based calculations 
are unrealistically large at ultra- high energies. 

TABLE IV: Charged current u-N cross sections, in cm'^, as a function of E^,, the laboratory energy, in GeV: ctbdhm from this 
work ; gctw from [s^l ; crcsMS from [s^ ; acQRS from 0|. Note that CTW and GQRS present values only up to = 10^^ 
GeV, while CSMS give values up to 5 x 10^^ GeV. The entry at Ei, = 10^^ GeV is marked by an asterisk to signal that the 
energy for asMS is actually 5 x lO^'^ GeV. 



E, (GeV) 


crsDHA/(cm^) 


acTW (cm^) 


ucsMs{cm^) 


CTGQflS (cm^) 


10*^ 


7.40 x 10"^" 


7.2 X lO"^'* 


6.9 X 10"^" 


6.34 X 10"^* 


10' 


2.07 x 10"^=^ 


2.0 X 10"^^ 


1.9 X 10-=^^ 


1.75 X 10"^^ 


10* 


4.89 X 10"^^ 


4.8 X 10"^^ 


4.8 X 10"^-' 


4.44 X 10"^^ 


10« 


1.00 X 10"^^ 


1.1 X 10"^^ 


1.1 X 10"^^ 


1.05 X 10"^^ 


10^" 


1.82 X 10"^^ 


2.2 X 10"^^ 


2.4 X 10"^^ 


2.38 X 10"^^ 


10" 


3.02 X 10"^^ 


4.3 X 10"^^ 


4.8 X 10"=*^ 


5.34 X 10"^^ 


10^2 


4.69 X 10"^^ 


8.3 X 10"^^ 


7.5 X 10"=*^* 


1.18 X 10"^^ 




14 



TABLE V: Neutral current v-N cross sections, in cm^, as a function of E,j, the laboratory energy, in GeV: (tbdhai, from this 
work; acTW, from [53|; crcsMS from [s^; acQRS from 0]. Note that CTW and GQRS only evaluate values up to Ei, — 10^^ 
GeV, while CSMS give values up to 5 x 10^^ GeV. The entry at = 10^^ GeV is marked by an asterisk to signal that the 
energy for asMs is actually 5 x 10^^ GeV. 



El, (GeV) 


o-BDHM (cm ) 


acTW [era ) 


acsMS[CTa ) 


o-GQRS (era ) 


lU 


O.U4 X iU 


o 1 n — 34 
Z. / X iO 


o c 1 n — 34 
2.0 X iO 


2. DO X iU 


10^ 


8.17 X 10"^'' 


7.6 X 10"^" 


7.3 X 10"^" 


7.48 X 10"^* 


10* 


1.97 X 10"^^ 


1.9 X 10"^^ 


1.9 X 10"^^ 


1.94 X 10"^^ 


10" 


4.12 X 10"^^ 


4.3 X 10"^^ 


4.4 X 10"=*^ 


4.64 X 10"^^ 


10i° 


7.58 X 10"^^ 


9.0 X 10"^^ 


9.6 X 10"^^ 


1.07 X 10"^^ 


10" 


1.27 X 10"^^ 


1.8 X 10"^^ 


2.0 X 10"^^ 


2.38 X 10"^^ 


10^2 


2.00 X 10"^^ 


3.5 X 10"^^ 


3.1 X 10"=^^* 


5.20 X 10"^^ 



C. Analytic form of the CC and NC cross sections, as a function of E^. 

In Fig. Owe plot the CC (upper points) and NC (lower points) cross section points of Table Hill cr, in rnb vs. E^, 
in GeV. We fit a 4 parameter fit, of the form cr = a + 61n_E,y + cln£'^ + din E'^, to the 12 points of Table Hill to obtain 
the analytic cross section functions 

(^cciK) = -19.91 + 4.685 Ini;^- 0.3798 ln2£;^ + 0.01078 In^i;^, (46) 
<yNciK) = -10.161 + 2.3041ni;^- 0.1801 ln2£;i. + 0.004926 In^-B^, (47) 

with Eu in GeV and the cross sections in mb. 




IQ6 IQlO JQl2 ^Ql4 JQl6 

Ev (GeV) 



FIG. 5: Plots of u-N cross sections, in cm'^, vs. E^, the laboratory neutrino energy, in GeV. Our analytic GC cross section 
is the upper solid curve of Eq. (|46p and the analytic NG cross section is the lower dashed curve of Eq. (|47[). The points are the 
numerical calculations of Table IIIII The agreement is better than 1 part in 1000. 

The upper solid curve in Fig. [5] is the In'^ E,^ parameterization of (Tqc md the lower dashed curve is the In^ E„ 
parameterization of (Tnc- The points are the numerical calculations of Table Hill Clearly, the In'^ E^, parameterization 
is excellent, with numerical agreement better than 1 part in 1000. 
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A full discussion of why a strong interaction Froissart bound of ln^(l/a;) gives rise to a weak interaction v — N 
cross section bound of In^ is given in the last paragraphs of Section IIIDI Conversely, a weak cross section u — N 
bound of In'^ E^, implies a strong cross section hadron-nucleon Froissart bound of In'^ s ~ In'^ W'^ , whose implications 
for hadronic physics we discuss next. 

D. Implications for hadronic physics 

We remark finally on the implications of these results for hadronic physics. Our fundamental assumption, discussed 
in detail in Part I [23, is that the virtual gauge boson-hadron scattering processes 7*p, W*N and Z*N are basically 
hadronic in nature, each interaction having the same Froissart-boundcd structure seen in all very high energy hadronic 
cross sections, including real 7p scattering (33-l34|. The structure functions F^^, F^'"', and FQ'^^'^^ for ep and CC 
and NC v{v)N scattering for a virtual boson mass = — are all related, differing only in the electromagnetic 
and weak charges and the helicity structure of the quark currents. Since, as we showed, F^^'^'^ is given to good 
approximation as a numerical multiple of F^^ with the multiplier a ratio of weak to electromagnetic charges, we 
expect the same Froissart-boundcd structure in neutrino as in electron interactions. However, the final hadronic cms 
energies W = potentially accessible and measurable in UHE neutrino interactions range far above the energies 
which have been studied in other experiments, 7 TeV at the Large Hadron Collider [28l - [30j and 57 TeV at the Pierre 
Auger cosmic ray array (3l| . This is shown in Fig. [6l where we plot the average value of the final cms hadronic 

energy W = = \/ /x versus the incident neutrino energy, with the average taken over the i'"^^''^ contribution 
to the CC neutrino cross section. It is evident from this log- log plot that Wav increases nearly as a power with E,y, 




IQfi 10* 10'" 10'2 lO''* lO"" 



Ey, GeV 

FIG. 6: Plot of the average cms hadronic energy W = \/S — \J jx in v-N scattering as a function of the incident neutrino 
energy. 

approximately as E^-'^^, close to the behavior W,nax ~ \/2mEy attained at the upper kinematic limit. 

The average cms hadronic energy in v-N scattering at E^ = 10^^ GeV, 112 TeV, is already twice the 57 TeV 
reached by the Auger Collaboration. As seen from Fig. |4] and Tables IIVI and |Vl the effects of the Froissart bound are 
clearly evident at this energy, with our predicted cross sections substantially lower than those obtained in standard 
approaches based on the DGLAP evolution of quark PDFs, extrapolated to small x. 

Experimentally demonstrating that v — N cross section measurements in this energy region are bounded by In'^ E^ , 
as shown in the preceding Section, would give a clear test of the entire picture discussed here. It would show that 
the virtual boson-hadron interactions are hadronic in nature, having the same Froissart-bound saturating structure as 
observed in other hadronic interactions, and thus allow future experiments to use weak v — N interactions to explore 
strong hadronic interactions at otherwise unobtainable energies. 
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IV. SUMMARY AND CONCLUSIONS 



In response to the ongoing need for theoretical calculations to guide design of ultra-high energy neutrino telescopes, 
we have improved our earlier calculations of the UHE v-N total cross section given in [s^ . The results include the 
b-quark contribution, whose importance we have shown in our analysis of the kinematical region that dominates the 
total cross section integral. The results are based on the all-orders relation between the dominant neutrino structure 
function F^*-'^-' and the 7p structure function established in Part I [l^l, and corresponding relations for the 
subdominant structure functions F!^''^^ and f^^'^-' correct to NLO. 

The basic assumption, justified in Part I, is that F^^, the 7*p reduced cross section for the interaction of an off-shell 
photon with a proton, has the Froissart-bounded form found in hadronic scattering and 7p scattering, which provides 
an excellent fit to the HERA data on ep DIS. The form can be extrapolated reliably through the factor ~ 2.5 extension 
of the range of the natural variable ln(l/a;) needed to reach the x values appropriate to UHE neutrino scattering. 

The uncertainties in our results arise primarily from the uncertainties in the values of the parameters in our saturated 
Froissart bounded fit to F^^^x, Q^) using the combined HERA data [5l|; they amount to 1% - 2% uncertainty in our 
cross section values over the entire E,j range that we study. We show that our v — N cross sections are bounded by 
In^ E^, which is a consequence of our Froissart-bounded In"" fit to ^^^(x, Q^). 

Comparing our cross section values to those of the most recent PDF-based study [s^ — which also includes the 
effects of the b-quark — we find that our cross sections are significantly lower than theirs at high neutrino energies, 
with the results diverging strongly for neutrino energies above ^ 10"'^^ GeV. 

In conclusion, we believe the neutrino cross sections calculated starting from our saturated Froissart bound fit 
to the existing experimental HERA data are the most physically motivated, and thus provide the best estimates of 
UHE energy v-N cross sections now possible. Moreover, these UHE neutrino interactions have significant potential 
implications for hadronic physics up to an average ultra-high hadronic cms energy = VI > 70, 000 TeV, if we can 
discover and measure v — N cross sections with laboratory energies E^, > 10^^ GcV. 
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